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We develop a coarse-graining procedure for two-dimensional models of fluctuating loops 
by mapping them to interface models. The result is an effective field theory for the scaling 
limit of loop models, which is found to be a Liouville theory with imaginary couplings. 
This field theory is completely specified by geometry and conformal invariance alone, 
and it leads to exact results for the critical exponents and the conformal charge of loop 
models. A physical interpretation of the Dotscnko-Fateev screening charge is found. 

1. Introduction 

Loop models are defined by drawing closed loops (which can come in one or more 
flavors) along the bonds of a two-dimensional lattice £, subject to the constraint 
that only loops of different flavor may cross. Boltzmann weights of different loop 
configurations are completely determined by specifying the weight (p) of all possible 
loop arrangements at a single vertex of C, and the fugacity (n) of every loop flavor. 
The partition function can be written as: § 

Z = Y,PT 1 PT---Pv V ^n l i...n l - , (1) 
e 

where F is the number of loop flavors and V the number of allowed vertex con- 
figurations. The sum in Eq. ([!]) goes over all allowed loop configurations Q; rrii 
(i = 1, . . . , V) is the number of vertices of type i, and lj (j = 1, . . . , F) the number 
of loops of flavor j in Q. 

The motivation for studying loop models comes from the observation that they 
are rather ubiquitous in the realm of two-dimensional statistical mechanics. Loop 
models appear very naturally when one considers domain boundaries in discrete spin 
models. For example, model A in Fig.l is equivalent to the Q = n 2 Potts model 
on the square lattice, where the loops run along the bonds of the medial lattice so 
as to encircle clusters of equal Potts spin. □ In a somewhat different setting, loop 
models appear as space-time diagrams of certain one-dimensional quantum systems, 
where the quantum partition function can be written as a sum over loop configu- 
rations with appropriate weights. This loop representation of quantum models has 
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been successfully taken advantage of in recent Monte Carlo (MC) simulations of 
spin chains and ladders .0 The loop representation allows for non-local (loop) up- 
date moves in the MC algorithm, which lead to a considerable reduction in critical 
slowing down, and rather large system sizes can be simulated. On the analytic 
side, the loop representation of the Heisenberg spin chain has lead to a very nice 
geometrical interpretation of correlation functions!] For instance, the spin-spin cor- 
relation function is associated with the probability that two points on the lattice 
are connected by a loop segment, in the appropriate loop model. 
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Fig. 1. Allowed vertex configurations in fully packed loop models on the square and honeycomb 
lattice; every vertex is assigned the weight p = 1. The models A and B have a single loop flavor 
with fugacity n, while loops in model C come in two flavors (black and grey), with fugacities n\ 
and U2- 

For a specific choice of Boltzmann weights, loop models can be in a critical state 
characterized by power law correlations. Here we examine the critical properties of 
fully packed loop models. These loop models have non- vanishing and identical (pi = 
1) vertex weights only for configurations that satisfy the fully-packing constraint; 
this constraint ensures that every vertex of the lattice is covered by one, and only 
one loop of each flavor. We have studied three such models, shown in Fig.l. Model 
A (No. of loop flavors, F — 1) is equivalent to the loop decomposition of the critical 
Q = n 2 Potts model,cl model B (F = 1) is the zero-temperature phase of the 
0{n) model, l3'tr0 and model C (F — 2) has recently been introduced as the loop 
generalization of the four-coloring model on the square lattice. □ The theoretical 
challenge we are faced with, is to find an effective description of a critical loop 
model in terms of a (conformal) field theory, from which critical exponents and 
other universal quantities can be calculated exactly. 

The answer to the above stated theoretical challenge is provided by mapping 
a loop model to an interface model, where loops are identified with contour lines 
of the interface. Coarse-graining of the interface model leads to an effective field 
theory which describes the critical fluctuations of loops in terms of a fluctuating 
height field. This effective field theory is equivalent by a duality to a Coulomb gas of 
electric and magnetic charges, with an additional electric charge placed at infinity. 

In the Coulomb gas approach to two-dimensional critical phenomena, the critical 
phase is identified with the vacuum phase of a gas of electric and magnetic charges;u 
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the strength of the Coulomb interaction is parameterized by the coupling constant 
g. Critical exponents are readily calculated once the value of this coupling is known. 
In the interface representation of loop models, g has the physical interpretation of 
the effective surface tension, or stiffness (see Eq. The effective stiffness is not a 
universal quantity and to calculate it in terms of the microscopic weights is typically 
as hard as solving the model exactly." In this paper we show that the stiffness of 
the interface is completely specified by geometry and conformal invariance alone. 
More precisely, it is shown that demanding conformal invariance of the effective 
field theory of critical loop models completely fixes all the couplings in the theory. 

2. Interface representations of loop models 

The general procedure for constructing an effective field theory of a loop model is 
outlined here for tha A model. Details of this construction for models B and C can 
be found elsewhere !M3 

Model A is the well known polygon representation of the Q = n 2 Potts model, 
and it is critical for < n < 2% For n > 2 the model has a finite correlation length, 
roughly the size of the largest loop; in the limit n — > oo it orders in one of two 
states, which are related by a lattice translation, and in which all the loops are of 
length four. 

The construction of an effective field theory for the scaling limit of a loop model 
can be broken up into three steps: 

2.1. Orient the loops 

Each loop is assigned an orientation by placing arrows along the bonds of the square 
lattice. In order to assign proper weights to the loops (n < 21 every left (right) turn 
is weighted by the phase factor exp(— ien7r/4) (exp(ien7r/4)).El For every closed loop 
on the square lattice the difference between the number of right, and the number 
of left turns is ±4. Therefore, by summing over the two possible orientations each 
loop is assigned the weight 

n = 2cos(7reo) . (2) 

2.2. Map to an interface model 

Each loop is viewed as a contour line of an interface model. For the A model 
this is the well known body-centered solid-on-solid (BCSOS) model;L3 for the B 
and C modeliLjnterface models with two and three component height variables 
are obtained.E3lrlL3 The microscopic heights z, of model A, are defined on the dual 
lattice, and z increases (decreases) by 1/2 every time an oriented loop is crossed 
from left to right (right to left). Coarse-graining of the microscopic height leads to 
an effective field theory of the loop model. 

a In the past, the value of the coupling constant has been determined for models that map to exactly 
solvable models by comparing the Coulomb-gas expression for a particular exeonent, which is g 
dependent, with the value of the same exponent found from the exact solutionjj 
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2.3. Coarse grain the interface model 

The interface model is coarse-grained in terms of the ideal states. These states 
minimize the variance of the microscopic height, i.e., they are macroscopically flat. 
In the loop model these are states with the maximum number of loops on the 
lattice, or equivalently, the states which one finds in the limit n — > oo. For finite 
n the lattice breaks up into domains of ideal states, and each domain is assigned 
a coarse-grained height equal to the average microscopic height over the domain: 
h =< z >. The continuum limit is taken by replacing the discrete heights by a 
height field h(x). The height field h(x) is compactified to a circle, i.e., 

h=h+K, K = {0,±1,±2,...} (3) 

where 1Z is the repeat lattice; the vectors in this lattice are height differences between 
equivalent loop configurations. 

The Euclidean action of the effective field theory, which describes the long-wave- 
length fluctuations of the microscopic height (or the critical fluctuations of the loops 
in the A model), can be written in terms of the field </>(x) = 2%h(x.) as 

Sg = 

Sb 

Sp = 

The three parts to the action 

(i) S g describes the Gaussian fluctuations of the height around the flat ideal 
states. 

(ii) Sb is the coupling of the height to the scalar curvature R. Namely, if a loop 
winds around a point of non-zero curvature, then the difference between the 
number of left and the number of right turns, for this loop, is no longer four; 
Sb corrects for this. For the square lattice in the infinite plane this term 
corresponds to a background charge 2eo placed at the point at infinity. 

(iii) S p is the potential term. Its origins are twofold: 1) it accounts for the discrete 
nature of the heights, and 2) it assigns proper weights (n < 2) to the loops, i.e., 
the operator w(<fr) is conjugate to the loop weight. From Eq.(^) we conclude 
that w{4>) can be expanded into a Fourier series, 

W{4>) =Y, W E ^ > ^ 
EeTZ* 

where 1Z* is the lattice dual to 1Z; here 1Z* is simply the set of integers. b Only 
the most relevant vertex operator (i.e., exponential of the height field) appear- 

b Elements of TV are the electric charges, while elements of Tt are the magnetic charges of the 
two-dimensional Coulomb gas, which is dual to the interface model. 



= Sgfr] + S b [</>] + S p [<l>] 

9 'dSc[V<^ 



= ^ d 2 xR, 



4tt 
iep 

ld 2 xw((t>) (4) 
have a simple geometrical interpretation: 
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ing in the above Fourier expansion is kept in the action. Upon examining 
the values which w(<j>) takes in the ideal states, we find this to be the vertex 
operator with electric charge 

E w = 2 . (6) 



3. Calculation of the coupling 

In order to calculate the value of the coupling g, for model A in the critical phase 
(0 < n < 2), we demand that the action S describe a conformally invariant field 
theory. This, as will be made clear below, is equivalent to the assumption that the 
loop fugacity does not change with the change of scale. 

For S to describe a conformal field theory the potentia l te rm has to be marginal, 
i.e., the scaling dimension x w of the operator w((f>) is two.ta The dimensionx^, can 
be expressed in terms of the coupling g and the background charge eo as:H3 

x w =E w (E w -2e )/2g = 2. (7) 

Using Eq. ([j]) we find the exact value of the coupling 

9 = 1 - e ; (8) 

both the background charge and g are continuous functions of the fugacity n, Eq.(||). 
This value of the coupling agrees with the value found from the exact solution of 
the eight-vertex model.H We emphasize that eo(n) in Eq.(||) defines a whole family 
of critical models, whose conformal charges are given in Eq.(^|). 

4. Discussion 

The field theory described by the Euclidean action in Eq. (||) is a Liouville theory 
with imaginary couplings, which has been suggested as the Lagrangean description 
of a two-dimensional Coulomb gas in the presence of a background charge.til The 
potential term S p is the screening charge, which was originally introduced into the 
Coulomb gas to ensure that the four-point correlation functions do not vanishJlj 
Here this operator appears quite naturally in the action as a result of the coarse 
graining, and it is the operator conjugate to the loop fugacity. 

The term S p also contains the locking potential, which enforces the condition that 
the heights are discrete. We find that the locking potential is marginal along the 
whole critical line (n < 2) of the loop model. This leads to the somewhat surprising 
conclusion that the associated interface model is at the roughening transition for 
all n < 2, not just at the boundary at n = 2. In the case of model B Baxter 
arrived at the same conclusion from the exact solution of the related three-coloring 
model; he showed that the partition function has a line of essential singularities for 
< n < 2.00 

For the B and C models the height field has more then one component and con- 
sequently the magnetic (M) and electric charges (E) are lattice vectors in 1Z, and the 
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dual 1Z* . The proper identification of these lattices follows from the coarse-graining 
procedure. The fact that the height is compactified on 1Z, Eq.(^), leads to some 
interesting conclusions about the symmetry of these loop models in the continuum. 
Namely, for n = 2 (B model) and n\ = = 2 (C model) the background charge 
vanishes, and these loop models are described, in the continuum, by the SU (3)fc = i 
and the 5't/(4)/ c= i Wess-Zumino-Witten model respectively:li-DE2l for the A model, 
at n = 2 one finds the SU(2)k=i WZW modelEj'Ej The action written in terms of 
the height field is the free-field representation of these sigma models. 

The conformal charge and all the critical exponents of the loop models can be 
expressed in terms of the coupling g and the background charge eo:El 

c = k 

9 

x(E,M) = J-E-(E-2e ) + f M 2 , (9) 
2g 2 

where x(E, M) is the scaling dimension of the operator with electro-magnetic charge 
(E,M), and k = 1,2, and 3 for the A,B, and C model respectively. For example, 
the m-RSOS models which are described in the continuum by the minimal models 
of conformal field theory can be mapped to the the A model with eo = l/(m + 
1), (m > 2).@ The conformal charge that we calculate from Eqs. (||) and (|J) is 
c = 1 — 6/m(m + 1), which is the well known expression. 
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